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Abstract 

By using Jackson's q-exponential function we introduce the generating 
function, the recursive formulas and the second order q- differential equa- 

r/\ , tion for the q-Hermite polynomials. This allows us to solve the q-heat 

equation in terms of q-Kampe de Feriet polynomials with arbitrary N 
moving zeroes, and to find operator solution for the Initial Value Problem 
for the q-heat equation. By the q-analog of the Cole-Hopf transformation 
we construct the q-Burgers type nonlinear heat equation with quadratic 
dispersion and the cubic nonlinearity. In g — >■ 1 limit it reduces to the 
standard Burgers equation. Exact solutions for the q-Burgers equation in 
^ ' the form of moving poles, singular and regular q-shock soliton solutions 

CO , are found. 






in 

^. ! 1 Introduction 

in ^ 

^D \ It is well known that the Burgers' equation in one dimension could be linearized 

^^ ■ by the Cole-Hopf transformation in terms of the linear heat equation. It allows 

one to solve the initial value problem for the Burgers equation and to get exact 
solutions in the form of shock solitons and describe their scattering. In the 
present paper we study the q-difFerential Burgers type equation with quadratic 
dispersion and the cubic nonlinearity, and find its linearization in terms of the 
C^ \ q-heat equation. In terms of the Jackson's q-exponential function we introduce 

the q-Hermite and q-Kampe-de Feriet polynomials, representing moving poles 
solution for the q-Burgers equation. Then we derive the operator solution of 
the initial value problem for the q-Burgers equation in terms of the IVP for the 
q-heat equation. We find solutions of our q-Burgers type equation in the form 
of singular and regular q-shock solitons. It turns out that static q-shock soliton 
solution shows remarkable self-similarity property in space coordinate x. 
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2 q-Exponential Function 

The q- number corresponding to the ordinary number n is defined as, [1], 

r 1 g"-l 



(1) 



where g is a parameter, so that n is the hmit of [n]q as q — > 1. A few examples 
of q-numbers are given here: [0], = 0, [1], = 1, [2], = 1 + q, [3], = 1 + q + q^. 
In terms of these q-numbers, the Jackson q-exponential function is defined as 

■A a;" 

For g > 1 it is entire function of x and when q — )■ 1 it reduces to the standard 
exponential function e^. The q-exponential function can also be expressed in 
terms of infinite product 



°° 1 1 

^^^"^ = n (1 - (1 _ q)qn^) = (l-(l-gX' (^) 

when q < 1 and 

e.(^) = n fi+(i--)4r^) = (i+i^-->Y , (4) 



q q / \ q 
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when q > \. Thus, the q-exponential function for q < 1 has infinite set of poles 
at 

Xn = ; r, n = 0, 1,.. (5) 

and for q > 1 the infinite set of zeros at 

a^« = - . _^y n = 0, 1,.. (6) 

The q-derivative is defined as 

Dqf(x) = — , (7) 

and when g — >■ 1 it reduces to the standard derivative Dqf{x) -^ J'{x). Using 
the definition of the q-derivative one can easily see that 

D-^{ax^) = a[n\qX^-\ (8) 

D^eq{ax) — aeq{ax). (9) 



3 q-Hermite Polynomials 

We define the q-Hermite polynomials according to the generating function 

^-^ j-n 

e,i-t^)e,i[2],tx) ^Y.H„ix;q)y—. (10) 

From this generating function we have the special values 

H2„(0;g) = (-ir^, (11) 

H2n+iiO;q) = 0, (12) 

where [n\ql = [l]q[2]q...[n\q, and the parity relation 

Hni-x;q) = {-iy'Hnix;q). (13) 

By q-difTerentiating the generating function pUj) according to x and t we have 
the recurrence relations correspondingly 

D.,Hn{x:q) = [2]q[n]qH,,^i{x-q), (14) 

i/„+i(x;q) = [2\qxHn(x;q) - [n]q Hn-i{qx;q) - [n]qq^~ Hn-i{y/qx;q). (15) 
Using operator 

Mq^q-^, (16) 

so that 

Mqfix) = f{qx), (17) 

relation (|15l) can be rewritten as 

Hn+l{x;q) - [2]qxH^{x-q) - [n],(A//, + g^Af^)iJ„_i(x;q). (18) 

Substituting ^ to dUl) we get 



H^+i{x-q) = \[2]qX~ ^ ^i?J i/„(a;;g) (19) 

By the recursion, starting from n = and i/g (2^) = 1 we have next representation 
for the q-Hermite polynomials 

A / Mq + q'^M fq \ 

Hn{x;q) = n f[2],^- ' ^^j^ -^ dA . 1 (20) 

We notice that the generating function and the form of our q-Hermite polyno- 
mials are different from the known ones in the literature, [5], [J], [3], [5]. 
In the above expression the operator 

Mq + q^ M^ = 2q'i qi''^ cosh[{lnqi){x- n)] (21) 



is expressible in terms of the q-spherical means as 

COsh[{\nq)x^]f{x) = i(/(qx) + f{-x)). (22) 

Using definition, [T], 

{x — ay^ — {x — a){x — qa) ■ ■ ■ {x — (7"^ a), 71 = f , 2, .. 

whicli now we apply for operators, we should distinguish the direction of multi- 
plication. We consider two cases 

{x - a);'< = {x- a){x -qa)---{x- q^'^^a), (23) 

and 

{x - a);'> = ix- q'^'^a) ■ ■ ■ {x - qa){x - a). (24) 

Then, we can rewrite (|20|) shortly as 

MgD iM^D^Y 



HUx;q)^m,x-^-^^)-q^^^^ .1 

First few polynomials are 

Ho{x;q) = 1, 

Hi{x;q) = [2],a;, 

H2ix;q)^[2]lx^^[2]„ 

H3{x;q) = [2]lx'-[2]l[3],x, 

m{x-q) = [2]\x^-[2]l[iUA],x^ + [2],[3] J2],.. 
When (? — > 1 these polynomials reduce to the standard Hermite polynomials. 

3.1 q-Differential Equation 



Applying D^ to both sides of (|T9)) and using recurrence formula ([14)) we get 
q-differential equation for q-Hermite polynomials 

^-^D^(Mq + q'^M^)D^Hn{x;q) - [2]gqxD^Hn{x;q) + [2]q[n]gqHn{x: q) = 0. 



3.2 Operator Representation 

Proposition 1 We have next identity 

e, (- [^^') e?([2],a;t) = eg{-t^)eg{[2]gxt). (25) 

Proof 1 By q- differentiating the q-exponential function in x 

Dle,{[2],xt) = {[2]tre,{[2\,xtl (26) 

and combining then to the sum 



°° [2]^a"i2» 



E T-n^reM^-t) = Y: '-\r^ e,i[2],xt\ 



a 
we have relation 



(27) 



e,(ai?^)eg([2],xi) = e,([2],at^)e5([2],a;t). (28) 

iJy choosing a = —l/[2]q we get 



e, -t:^^' e,([2],xt) = eg{^t')eg{[2]gxt). (29) 



1 

[2], 

Proposition 2 T/ie neif identity is valid 

Hr.{x;q) = [2re,(^-^^Diy-. (30) 

Proof 2 r/ie rz(;/it /lanrf sirfe of li25\) is the generating function for q-Hermite 
polynomials. Hence expanding both sides in t we get the result. 



Proposition 3 



-■q 



[2],r "[2], l^^' [2], '\ [2 



Proof 3 we use i30\) and relation U9\) . 

Corrollary 1 If function f{x) is analytic and expandable to power series f{x) 
Sn=o '^"2;" then we have next q-Hermite series 



■-q 



^^i)H')-t''-^M§^^ (32) 
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4 q- Kampe-de Feriet Polynomials 

We define q-Kampe-de Feriet polynomials as 

H„{x,t;q) = {-iyt)^Hn(j-^^^=), (33) 

and from p9p we have the next recursion formula 

H„+i{x,t;q) = [x + {Mq + q^M^)vtD^^ Hn{x,t;q). 
By recursion it gives 

n 

Hn{x, t;q) = l[(^x+ (M, + q^M^^tD,^ ■ 1 (34) 

or by notation (IMl) 

Hn{x,t-q) = {{x + MqVtD^)+q^M^vtDX -1 

First few polynomials are 

HQ{x,t;q) = 1, 

Hi{x,t;q) = a;, 

H2{x,t;q) =x^ + [2]gj^t, 

H3{x,t;q) = x^ + [2]q[3]qiytx, 

Hi{x,t; q)^x^ + [3] J4], j.ia;2 + [2] j3]q[2],2;/2^2^ 

5 q-Heat Equation 

We consider the q-heat equation 

{Dt~iyDl)(l3{x,t)=0. (35) 

Solution of this equation expanded in terms of parameter k 

(l){x, t) = eq{vk'^t)eq{kx) = ^ — _ff„(x, i; q), (36) 

n=0 ^ -I' 



gives the set of q-Kampe-de Feriet polynomial solutions for the equation. Then 
we find time evolution of zeroes for these solutions in terms of zeroes Zk{n, q) of 
q-Hermite polynomials, 

Hn{zk{n,q),q)^0, (37) 

so that 

Xk{t) = [2]zk{n,q)^^Ui. (38) 

For n=2 we have two zeros determined by q-numbers, 



a;i(i) = ./[2],V=j4, (39) 



a;2(i) = -^[2]gV^l4, (40) 

and moving in opposite directions according to (I38p . For n=3 we have zeros 
determined by q-numbers, 



xi{t) = -^[3]glV^^t, (41) 

X2{t) = 0, (42) 

xsit) = yjsy^/^^, (43) 

two of which are moving in opposite direction according to psp and one is in 
the rest. 



6 Evolution Operator 

Following similar calculations as in Proposition I we have next relation 

Cq {vtDl) eq{kx) = eq{vtk^)eq(kx). (44) 

The right hand side of this expression is the plane wave type solution of the 
q-heat equation 

{Dt-i^Dl)4>{x,t)=Q. (45) 

Expanding both sides in power series in k we get q-Kampe de Ferie polynomial 
solutions of this equation 

Hn{x,t;q)=eg{i^tDl)x'\ (46) 

Consider an arbitrary analytic function f{x) — X]n=o'^"'^"' then function 



f{x,t)^eq{vtDl)f{x) = ^ a„e, (^.ii^l) x" (47) 

oo 

= ^a„i7„(x,t;<7), (48) 



n=0 



is a time dependent solution of the q-heat equation dJS]). 

According to this we have the evolution operator for the q-heat equation as 

U{t) = eg [vtDl) . (49) 

It allows us to solve the initial value problem 

{Dt~vDl)dp{x,t) = 0, (50) 

0(x,O+) = fix), (51) 

in the form 

0(x, t) = e, {vtDl) 4>{x, 0+) = e, {vtDl) f{x). (52) 

7 q-Burgers' Type Equation 

Let us consider the q-Cole-Hopf transformation 

u{x,t) = -2v , (53) 

(/)(a;,i) 

then M(a;, t) satisfies the q-Burgers' type Equation with quadratic dispersion and 
cubic nonlinearity 

Dtu{x, t) - Dlu{x, t) = i [{u{x, qt) - u{x, t)M^)D^u{x, t)] ~ 

i [Dx {u{qx, t)u{x, t))] + J \u{q^x, t) — u{x, gi)] u{qx, t)u{x, t). 
When (jf — > 1 it reduces to the standars Burgers' Equation. 

7.1 I.V.P. for q-Burgers' Type Equation 

Substituting the operator solution ([5^ to ([5^ we find operator solution for the 
q-Burgers type equation in the form 

This solution corresponds to the initial function 

uix,0+) = ~2:.^^. (55) 

fix) 

Thus, for arbitrary initial value problem for the q-Burgers equation with u{x, 0^) : 
F{x) we need to solve the initial value problem for the q-heat equation with ini- 
tial function f{x) satisfying the first order q-differential equation 

(D, + l.F{x))f{x) = 0. (56) 



8 q-Shock soliton solution 

As a solution of q-heat equation we choose first 

(p{x, t) = eq (fc i) Cq [kx) , 

tlien we find solution of the q-Burgers equation as a constant 

u{x,t) — —2uk. 
If we choose 

4}{x, t) = Bq [kit) eq (kix) + Eq [kjt) Cq {k2x) , 

then we have the q-Shock soliton solution 

klCq (kit) Cq (fcix) + k2eq (A:|t) Cq {k2x) 



u(x, t) — —2v- 



eq (kjt) Cq (kix) + eq [kit) eq {k2x) 



(57) 

(58) 
(59) 

(60) 



Due to zeroes of the q-exponential function this expression admits singularities 
for some values of parameters fci and fc2. In Fig.l we plot the singular q-shock 
soliton for fci = 1 and /c2 = 10 at time i = 0. 
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Figure 1: Singular q-Shock Soliton 

However for some specific values of the parameters we found the regular 
q-shock soliton solution. We introduce the q-hyperbolic function 



coshg(a;) 



eq{x) -f eq{-x) 



(61) 



or 



cosh,(x) = - e,(x) + — -- , (62) 

2 \^ e^^{x)) 

then by using infinite product representation for q-exponential function we liave 



coslig(x) = - 




From ([SI),® we find that zeroes of the first product are located on negative axis 
X, while for the second product on the positive axis x. Therefore the function 
has no zeros for real x and coshq(O) = 1. 

If fci = 1, and ^2 = ^Ij the time dependent factors in nominator and the 
denominator of ([BUI cancel each other and we have the stationary shock soliton 

M(x,t) = -2iy^44^^4^--2j^tanh„(x). (63) 

eq{x)+eq{-x) 

Due to above consideration this function has no singularity on real axis and we 
have regular q-shock soliton. 

In Fig. 2, Fig. 3 and Fig. 4 we plot the regular q-shock soliton for ki — 1 and 
fc2 = — 1 at different ranges of a;. It is remarkable fact that the structure of our 
shock soliton shows self-similarity property in space coordinate x. Indeed at the 
ranges of parameter x = 50, 5000, 500000 the structure of shock looks almost 
the same. 

For the set of arbitrary numbers fci, ..., kn 

N 
0(2;, t)^^e.q {kit) Cq {knX) , (64) 

n=l 

we have multi-shock solution in the form 

u{x, t) = _2^ ^^^i^"^^(^"^)^^(^"^) . (65) 

In general this solution admits several singularities. To have regular multi- 
shock solution we can consider the even number of terms N — 2k with opposite 
wave numbers. When A^ = 4 and ki = 1, fc2 = — 1,A:3 = 10,fc4 = —10 we have 
q-multi-shock soliton solution, 

eq{t) sm\[x) + 10e,(100t) sinh,(10x) 
"^^' ^ *" e,(t)cosh,(a;)-fe,(100i)cosh,(10x) ■ ^ ' 

In Fig. 5 we plot A^ = 4 case with values of the wave numbers fci = 1, 
A:2 = — 1,^3 = 10,A;4 = — 10 at i = 0. To have regular solution for any time t and 
given base q, we should choose proper numbers ki which are not in the form of 
power of q. This question is under the study now. 
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Figure 2: Regular q-Shock Soliton For fci — 1, k2 — —1, at range (-50, 50) 
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Figure 3: The regular q-shock soliton for ki = 1, fc2 = —1 at range (-5000, 5000) 
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Figure 4: The regular q-shock soliton for fci — 1, k2 = —1 at range (-500000, 
500000) 
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Figure 5: q-shock regular 
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